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Spatial gaps correspond to the projection in position space of the gaps of a periodic structure 
whose envelope varies spatially. They can be easily generated in cold atomic physics using finite- 
size optical lattice, and provide a new kind of tunnel barriers which can be used as a versatile tool 
for quantum devices. We present in detail different theoretical methods to quantitatively describe 
these systems, and show how they can be used to realize in one dimension matter wave Fabry- 
Perot cavities. We also provide experimental and numerical results that demonstrate the interest of 
spatial gaps structures for phase space engineering. We then generalize the concept of spatial gaps 
in two dimensions and show that this enables to design multiply connected cavities which generate 
a quantum dot structure for atoms or allow to construct curved wave guides for matter waves. At 
last, we demonstrate that modulating in time the amplitude of the periodic structure offers a wide 
variety of possible atom manipulations including the control of the scattering of an incoming wave 
packet, the loading of cavities delimited by spatial gaps, their coupling by multiphonon processes 
or the realization of a tunable source of atoms. This large range of possibilities offered by space 
and time engineering of optical lattices demonstrates the flexibility of such band gap structures for 
matter wave control, quantum simulators and atomtronics. 

PACS numbers: 03.75.Kk,03.75.Lm,67.85.-d 


I. INTRODUCTION 

Metamaterials are a fast-developing research field with 
a wide variety of applications [l]. A famous example in 
optics is provided by photonic band gap structures that 
exploit multi-wave interferences [2]. Optical components 
such as wave guides, bending light, micro-resonators or 
filters can be realized by introducing periodicity defects 
at suitably selected spots within a crystal. By pairing 
them it is now possible to design photonic circuits. 

Inspired by those developments in optics, a few stud¬ 
ies haveproposed to investigate their counterpart in atom 
optics p-0. The first element that was envisioned was 
the Bragg or multilayer dielectric mirror in which the 
layers were provided by a laser light interference pattern 
i.e. a periodic succession of dark and bright layers with a 
typical size of hundreds of nm II- In contrast with op¬ 
tics, the multiple matter wave interferences are not used 
here to improve the reflectivity but to gain on velocity se¬ 
lectivity and to design ultra narrow and tunable velocity 
filters [i, [6]. Indeed, many quantum devices require spe¬ 
cific shaping or momentum-selective filter of an atomic 
wave packet Ik!- 

The experimental realization of such a mirror was per¬ 
formed by studying the scattering of a guided Bose- 
Einstein condensate on a one-dimensional and finite- 
size optical lattice [18]: the band structure was directly 
probed and a wide variety of velocity filters (notch, band 
pass, high and low pass) has been demonstrated. The en¬ 
velope of the optical lattice projects in real space the gaps 


of the band structure. The corresponding spatial gaps 
generate a new type of atomic tunnel barriers. Cross¬ 
ing such a barrier in real space amounts to performing 
a Landau-Zener transition between two adjacent bands 
[H Eo|. In (2l|, a matter wave cavity delimited by 
two spatial gaps has been investigated. Among many 
prospects, the shaping of the envelope using for instance 
Spatial Light Modulators (see e.g. |22|) should offer the 
possibility to realize a mode-locked atom laser. 

In contrast with structured materials used in optics, 
the optical lattice depth can be modified in time. This 
possibility was exploited experimentally in (23[ to realize 
tunable velocity filters and mirrors. This brings forward 
a new tool for the emerging field of atomtronics US 
|4]: the components that will be combined to realize an 
atomic circuit could have a tunable functionality. 

In this article, we detail different methods to calcu¬ 
late in practice the properties of spatial gaps and pro¬ 
pose several possible applications of this concept, from 
the implementation and observation of complex quan¬ 
tum phenomena to the construction of versatile tools for 
atomtronics. To this aim, we study in detail different 
setups in one or two dimensions using static or time- 
modulated potentials. After a short reminder about the 
Bloch band formalism (Sec. EK) , the Hill’s method is pre¬ 
sented in Sec. Efc. This mathematical trick provides the 
exact imaginary part of a wave vector n inside the gaps 
of a one-dimensional, infinite and uniform periodic po¬ 
tential. The characterization of those evanescent matter 
wave modes plays a key role in the understanding of spa¬ 
tial gaps. However, the exact mathematical treatment is 
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restricted to one-dimensional system. In Sec.HIlG. we use 
those exact results to validate a systematic perturbative 
approach to evanescent modes that can be readily gener¬ 
alized to higher dimensions. Those results about evanes¬ 
cent modes are then used to define precisely the concept 
of spatial gaps within a locally uniform potential depth 
approximation and discuss their properties (Sec. ED). In 
Sec. US we discuss the dynamics of a wave packet in the 
presence of two spatial gaps (Fabry-Perot like device) or 
more. We also demonstrate experimentally the interest 
of such setups for phase space engineering. In Sec. IVl we 
address the generalization of the previous ideas to two 
dimensions leading to the design of simply and multi¬ 
ply connected matter wave cavities and of curved wave 
guides. Section [VI] is devoted to the generalization to 
the time domain of the concept of spatial gaps. We show 
how multiple cavities can be designed in particular in 
the wing of the envelope of the optical lattice and how 
their absolute and relative population can be controlled 
by multiphonon processes. 

II. GAPS AND EVANESCENT MODES IN AN 

OPTICAL LATTICE 

A. Bloch formalism 


conditions. As such, we expect a quantification of the 
energies that correspond to the discrete band index n. 
Because of the periodicity of the Bloch functions u n j c , 
we can limit ourselves to the first Brillouin zone, i.e. 
k E [—&r/2, &r/ 2]. By expanding the periodic Bloch 
functions u n as Fourier series, we get 

V’ra(*) = e ikz u n , k (z) = v k+ e kR e^ k+ik ^ z . (4) 

i 

For a given fc, working out the Bloch state i/j n ,k(z) 
amounts to finding the coefficients Vk+£k n • Similarly, the 
potential U can be Fourier expanded: 

U(z) = J2u P e ipkRZ . (5) 

P 

Combining Eqs. and © with the Schrodinger equa¬ 
tion ©, we infer the following set of coupled equations: 

h 2 

— (k + £k R ) 2 Vk+£k R + Up v k+(£-p)k R = Evk+ek K - 

p 

(6) 

To find the eigenvalues E n (k ), we truncate this infinite 
linear system and solve the resulting finite size matrix 
equation. Consider an attractive lattice of potential 


Before tackling the detailed treatment of spatial gaps, 
we first provide hereafter a short reminder on the Bloch 
formalism for quantum particles evolving in a periodic 
structure. Consider a particle of mass m in a one¬ 
dimensional periodic potential of period d, U(z + d) = 
U(z). The stationary Schrodinger equation that de¬ 
scribes the dynamics in such a potential reads: 

W{z) ={E+ U{z)^j ${z) = E'ip(z). (1) 

The Bloch theorem states that the eigenstates of a pe¬ 
riodic Hamiltonian may be written as the product of a 
plane wave function with wave-vector k and a function 
u n,k{ z ) that has the same periodicity as that of the po¬ 
tential U [Hj: 

1pn,k(z) = e lkz Un tk (z) with Un :k (z + d) = Un :k (z). (2) 

Except for k = 0, these Bloch states are propagating 
states. The eigenenergies associated with i/j n: k , E n {k ), 
are periodic with period k R = 2it/d\ E n (k) = E n (k+kn). 
By inserting W/c in Eq. ©, we deduce that the func¬ 
tions u U: k are the eigenstates of an effective /^-dependent 
Hamiltonian: 


C p + hk ) 2 

Hk 2m +U{Z) • (3) 

Using this eigenvalue equation and the boundary condi¬ 
tion u U: k(z) = u n ,k(z + d), the signification of the band 
index becomes clear. E n (k) are the eigenvalues of an 
effective Hamiltonian in a box with periodic boundary 


U(z) = -U 0 {cos{k R z) + l)/2, (7) 

where the Uo > 0 is the depth of the lattice. The matrix 
M that provides the eigenvalues through M • V = EV 
has a simple band structure with 


/ b-N u 

u b-jy+l u 


\ 


M = 




u 6/v-i u 

u bjsr ) 


( 8 ) 


where u = — Uo/4, bi = E R (k/k R + E) 2 — Uq/2 and 
Er = T?k\l2m. In this way, one finds for each k in 
the Brillouin zone, a discrete set of energies E n (k). By 
plotting these eigenvalues for all values of k , we obtain 
the energy bands that are separated by energy gaps. 
The border of the gaps are obtained either at k = 0 
or k = ±k R /2. In the energy gap there exist solutions 
of the eigen-equation for complex values of the wave vec¬ 
tor k. They are associated with evanescent modes and 
are usually invoked in solid state physics only when deal¬ 
ing with surfaces and junctions. They are at the heart 
of this article. In the next subsection, we provide a re¬ 
minder on Hill’s method to calculate the expression of 
those complex wave vectors. We then propose an alterna¬ 
tive method to calculate them approximatively through 
a perturbative approach. This latter method has the ad¬ 
vantage of being more physical and can be directly gen¬ 
eralized to higher dimensions. 
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B. Hill’s method 


Using the potential the Schrodinger equation is 
nothing but a Mathieu equation. Calculating the com¬ 
plex wave vector associated with an energy lying in a 
band gap amounts to calculate the so-called Mathieu 
characteristic exponent ft G C. In order to determine 
them, we detail hereafter the Hill method following [ 23 - 

H3. 

The eigenvalue matrix equation can be recast in the 
form A (ft; E, Uq) • V = 0 with 


A{k;E,U 0 ) 


1 t-N \ 

C-iV+l 1 £-iV+l 


6v-1 1 £/V-l 

V 6v 1 / 


(9) 


where & = (-U 0 /4£ R )/((ft/& R +f) 2 -C) and C = (U 0 / 2+ 
F^/Fr. Nontrivial solutions are obtained for A (ft) = 
det(A) = 0. In the limit N —>> oo, A (ft) is periodic with 
period &r. We can thus restrict the domain of Re [ft] to 
the domain: 0 < Re [ft] < &r. Since each £g(ft) appears 
only once per line, A is a sum of products of the £g(ft). 
The & are analytical functions of ft except at poles which 
are of finite order (meromorphic functions); A is thus 
analytical except at ft/£;R = a/C — Since each function 
appears only once in A, all the poles of A(ft) are simple 
and located at ft/fc R = y/Q — A 

The trick of Hill’s solution consists in considering the 
function [2j| 


cos(27Tft/A:R) — cos(27Tv/C) 

that has exactly the same poles as A. One can thus 
choose a function C(K) such that the function defined by 

0(ft) = A (ft) - C(ft) • £>(ft) (11) 

has no singularities. In the interval [0, fc R ], A has only 
one pole and the function C can be chosen constant and 
equal to the ratio between the residues of the functions 
A and D at this unique pole. With such a choice for 
C, the function O is analytical on the whole complex 
plane (holomorphic function). It must then be a constant 
according to Liouville’s theorem. 

To determine the value of ( 7 , we use the limit ft — »• 
+ioo. We find D(k —»• +ioo) = 0 and A(ft -> -Moo) = 1, 
and deduce that ©(ft— »+zoo) = l. As a result, we 
obtain the value of C = (A(ft) — 1 )/D(k). Using D{ 0) = 
1/(1 — cos(27rvD)? we infer II s explicit expression 

c = (A(0) - 1) • (1 - cos(2 TT^/C)). (12) 


From A (ft) = 0 (non-trivial solution), we finally get the 
equation fulfilled by ft for a given energy E (through the 
parameter (): 

cos(27rft/A:R) = 1 — A(0)(1 — cos(27r^/C))- (13) 


Real solutions of this equation provide the equations for 
the energy bands E(k ), while those with an imaginary 
part account for the evanescent modes in the gaps. 


C. The perturbative approach 

The previous approach is exact but requires to solve 
numerically the nonlinear equation m- It is possible 
to work out a perturbative treatment of the evanescent 
modes, for two important purposes: to get an explicit 
formula for the imaginary wave vector valid in the low 
lattice depth limit and to use it as a starting point for 
obtaining the imaginary wave vectors for 2D or 3D optical 
lattices. 

To investigate perturbatively the first gap, we use the 
two mode approximation for the potential ([7j). It consists 
in reducing the infinite set of equations to two coupled 
equations involving the modes at the edge of the Brillouin 
zone k = ±&r/2 J29[: 

(£r/4 - Uo/2)v kn / 2 - (Uo/E)v _ kR /2 = Evj , CR / 2 , 
-(Uo/4:)v kn /2 + (£r/4 - Uq/ 2)v_ kR /2 = Ev_ kK / 2 (14) 


Introducing the reduced depth parameter so = Uo/T^r, 
the gap borders deduced from the previous set of equa¬ 
tions are given by E± = (Ur/4)(1 — 2so ± so)- T° deter¬ 
mine the imaginary part of wave vectors inside the cor¬ 
responding band gap, we set the energy equal to E x = 
(Ur/4)(1 — 2so + so^) with — 1 < x < 1. We then search 
for a two-mode solution with k± = ±(1 =b iK){k^/2) 
where K depends on x and 5o: 

(E r (1 + iK) 2 - E r - xUo)vk K/2 - UoV-k K/2 = 0, 
-Uov kn / 2 + (^ r (1 - iK) 2 - E r - xUo)v_ kn / 2 =(®^) 


This set of equations has a solution if K obeys the equa¬ 
tion 


A 4 + K 2 {4 + 2 xs 0 ) + s 2 {x 2 - 1) = 0 (16) 

whose solution is given by 


K(x,so) = )J-2-s 0 x + ^/4 + s 2 +4s 0 x 

^ s 0 \/l - x 2 /2, for |s 0 | <C 1. (17) 

Figure [l] summarizes the comparison between the exact 
method and the perturbative approach detailed above for 
different lattice depths so = 0.1, so = 0-3 and so = 0.5. 
We observe that the maximum of the bell shaped curve 
for K is indeed on the order of so/2 as expected from 
the perturbative treatment. As expected the accuracy 
of the perturbative method based on the two mode ap¬ 
proximation decreases as so increases. In appendix A, we 
discuss more systematically the validity of the two mode 
approximation. 
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FIG. 1: Imaginary part of the wave vector normalized to 
/cr/2, K(E, so), in the first band gap for so = 0.1, so = 0.3 
and so = 0.5: exact calculation (solid line), perturbative cal¬ 
culation to the lowest order (dashed line) and perturbative 
calculation pushed to the next order (black square). 


The strategy to get a higher accuracy is quite clear, 
we have to add extra coupled modes. In the perturba¬ 
tive approach presented in the previous paragraph we 
have only taken into account the coupling between two 
modes v_ kK / 2 and v kn / 2 . To get a better approximation 
we can repeat the previous argument with four modes: 
^-3fc R /2, U-fca/2, v kn /2 and v 3kn/2 . The result is shown 
in Fig. [l] for so = 0.5. We already obtain a quite good 
agreement with the exact result (relative distance below 
1.7 %). Increasing further the number of coupled modes 
will increase the convergence towards the exact result. 
As expected for such a perturbative approach, to reach 
a given accuracy the number of modes that have to be 
taken into account increases with the potential depth (so 
parameter). We have also extended this approach to the 
second gap in Appendix B. 


a Landau-Zener transition between two adjacent bands 

0113 . 



FIG. 2: Bloch band distributions in the first Brillouin zone for 
two different potential depths s i and s 2 > s i. The pseudo¬ 
energy (dashed line) is fixed so that E+(si) = E-(s 2 ). 


To calculate the transmission probability T(E ), 
through the barrier, we shall calculate the value of the 
imaginary part Im[ft(z; E)] in the interval [zi;z 2 ]: 


T(E) = exp ^—2 J 


Im[tt(z; E)]dz 


(18) 


In this interval, the value of the depth s varies with space 
but the pseudo-energy, E, remains fixed. For s in the 
interval [si; s 2 ], the value of x is set by the energy: 


E hUP, 


where x depends on 2 through the local depth parameter 
s{z). At the scale of the width of the barrier, we can 
assume that the envelop varies linearly: s(z) = s(zi) + 
(z — zi) (d z s) . We deduce 


III. SPATIAL GAPS AS TUNNEL BARRIERS 
FOR MATTER WAVES 

The spatial variation of the envelope of an optical lat¬ 
tice projects the band gaps in real space. The latter are 
therefore referred to as spatial gaps. In the following we 
consider the first gap of a lattice for which the scale of 
variation of the envelope is large compared to the lattice 
period. Under this assumption, one can consider the lat¬ 
tice locally uniform at each position. For a given pseudo¬ 
energy E and potential depth s, we consider the closest 
gap at an energy below E. There are two values of the 
depth si < s and s 2 < s such that E+(si) = E-(s 2 ) = E 
with s\ = s(zi) and s 2 = s(z 2 ) (see Fig. [2j ). In the in¬ 
terval [z 1 ; z 2 \ , the Mathieu exponent k acquires an imag¬ 
inary part. The size of the corresponding spatial barrier 
is 8z = \zi — z 2 1. This virtual barrier thus acts as a reg¬ 
ular tunnel barrier in the energy range spanned by the 
gap. When crossing such a barrier, the atom undergoes 


where f(E) is a function that depends only on the energy: 

f(E) = 2 [ 2 Im[K(E,s)]ds. (21) 

J Si 

This exact form therefore exhibits a probability of 
transmission that decreases exponentially with the in¬ 
verse gradient of the envelope. The smaller the gradient, 
the larger the tunnel barrier width. This gradient there¬ 
fore appears as a parameter that can be tuned experi¬ 
mentally to set the transmission tunnel probability. 

This general formalism can be applied in combina¬ 
tion with the perturbative approach of Sec. Hit 1 to 
get explicit formulas for the transmission probabil¬ 
ity. In the perturbative limit, we have Im[ft(e, s)] = 
&r\A 2 (4e — 1 + 2s) 2 /4 (using Ea.flTTl)). si = (1 — 
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4e)/3 and S 2 = 3si (where e = E/Er), and find for 
e < 1/4. 


m = kR nvs {1 ~ ie) 


( 22 ) 


It is worth noticing that the energy dependence of this 
transmission probability is quite different from that ob¬ 
tained for a real repulsive potential barrier. This turns 
out to be an advantage of this kind of tunnel barrier 
since one can achieve pretty high transmission probabil¬ 
ity for a single barrier as experimentally demonstrated 


21j . This is a feature that is difficult, if not impossi¬ 


ble, to obtain using an optimally focussed blue detuned 
laser to realize a repulsive barrier for atoms because of 
the diffraction limit l30l. 


IV. SPATIAL GAP STRUCTURES IN ONE 
DIMENSION 

The preceding section considered a spatial gap as an 
isolated barrier. However, it is quite easy to engineer 
systems where several spatial gaps are present at different 
locations. In this section, we address the new features 
that emerge when one consider a pair of spatial gaps 
and discuss the case of multiple barriers. We also show 
that this type of device enables to engineer the phase 
space distribution of wave packets and demonstrate it 
experimentally. 


A. Matter wave Fabry-Perot cavity 

Consider a periodic potential with a Gaussian envelope 
as used experimentally in [18] (Fig. 03(a)): 

U(z) = - Vorsin 2 (^) (23) 

with Vq(z) = Uq exp (—2 z 2 /w 2 ). This Gaussian shape re¬ 
sults from the envelopes of the two Gaussian laser beams 
that generate the optical lattice. The parameter Uo > 0 
is proportional to the intensity of the laser beams and 
the lattice spacing d is proportional to the laser wave¬ 
length. Hereafter, we choose the typical experimental 
values w = 140 fi m, d = 0.65 /im and a potential depth 
Uo = F R /2. 

This potential is symmetric and therefore produces 
pairs of symmetric spatial gaps. In this section, we focus 
our study on a single pair of spatial gaps. Each spatial 
gap acts as a matter wave mirror of energy-dependent 
reflectivity (see Fig. [31(b)). As such, the potential expe¬ 
rienced by atoms is reminiscent of that of a Fabry-Perot 
cavity [2lj |. 

The complex amplitude, A, of the output wave is ob¬ 
tained by adding the contribution of the multiple reflec¬ 
tions of matter waves inside the cavity and reads (see 
Fig. ED a)) 

A = t 2 e iSvi + t 2 r 2 e iSV2 + t 2 r 4 e iS<fi3 + ... (24) 


with S(fi = 2 p t + p, Sip 2 = 2 pi + 3p + 2 p r , ... For sake 
of clarity, the letter t (resp. r) refers to the modulus of 
the transmission (reflection) amplitude, the correspond¬ 
ing phase is captured in the term p t (resp. p r ) [3lj |. 
To determine the transmittance function, Tpp = \A \ 2 , of 
this Fabry-Perot like device, we square the modulus of 
the amplitude: 


Tfp(F) = 


t 2 e 2i<pt e i<p(E) [ l + V^ 

V 

1 1 


1 + 4 


T 2 {E) T(E) 


r 2n e 2in(<p(E)+ip r ) 


sin 2 (<p(E) + (p r ) 


(25) 


The transmission probability Xpp(F) is related to the 
transmission coefficient by T(E) = t 2 (an explicit ex¬ 
pression is given in Eq. (jT8])) . The phase accumulated 
over the cavity length can be calculated by a semiclassi- 
cal approach: 

p(E) = ^ k dz = J R e[hv(z;E)\ dz. (26) 

We consider a single path from the first barrier to the sec¬ 
ond, this is the reason why a factor 1/2 appears in the 
previous formula. As expected, the phase of the trans¬ 
mission amplitude, p t , is washed out in the transmission 
probability. However, the phase of the reflection wave, 
p r , is responsible for a correction to the semiclassical 
formula that slightly depends on the energy (see below). 

In Fig. [31(c), we compare the transmission Xpp(F) ob¬ 
tained from our semiclassical approach with the trans¬ 
mission Xpw(F) obtained by solving numerically the sta¬ 
tionary Schrodinger equation using plane waves (referred 
to as the exact result in the following). To obtain an 
almost perfect agreement for the amplitude and without 
loss of generality, we have included a global 7r/4 offset 
in the phase of the semiclassical model. A more careful 
analysis reveals that a small energy drift of the position 
of the transmission peaks remains. The inset of Fig. He) 
represents the corresponding energy-dependent phase dif¬ 
ference, p r — 7r/4, between the semiclassical model and 
the exact solution of the Schrodinger equation (the shift 
correction is on the order of 1% (0.03 rad per peak)). 
This small residual shift confirms the good predictability 
of the semiclassical approximation. 

To perform a direct comparison of the energy depen¬ 
dence of the phase acquired in transmission through 
the cavity, we proceed in the following manner: (i) 
we identify precisely each transmission peak of the 
Schrodinger simulation and (ii) we assign a n phase 
shift between two successive peaks. In Fig. [H we plot 
the corresponding cumulated phase (black dots) and 
compare it with the semiclassical evaluation of the 
phase p(E ), performed through Eq. (l26| ) using the Hill’s 
method (see subsection EH). We find once again a 
very good agreement between the two methods [32l |. 
As intuitively expected, the phase strongly depends on 
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(a) 


-Z2-Z1 z 1 £2 



-200 -100 0 100 200 



E (E R ) 


FIG. 3: (Color online) (a) Optical lattice with a Gaus¬ 
sian envelope showing two spatial gap barriers in the spaces 
[—Z2,—zi\ and [ 21 , 22 ]. (b) Density plot of Im(fc) associated 
to the potential (1251 . as a function of the position and the 
energy. The solid black lines correspond to the values of the 
resonances of this cavity (see text), (c) Black line: Exact solu¬ 
tion for the transmission, Tpw, as a function of the energy E. 
Green/grey line: Transmission Tfp calculated by the semi- 
classical approach (I25jl including a global offset (ip r — 7t/4). 
Inset: Phase difference (ip r — 7r/4) between the semiclassical 
model and the exact one as a function of energy. 


the waist w of the gaussian envelope (Fig. (D^b)). It is 
therefore possible to tune the number of resonances by 
changing the cavity width, in accordance with what 
could be expected from standard phase space volume 


(a) 


T3 

03 

r—I 




(Er, 2 w) 




0.125 0.15 0.175 0.2 

E(E r ) 


0.225 0.25 


FIG. 4: (Color online) Solid green/grey line: numerical inte¬ 
gration of the semiclassical phase (l26l) . Black dots: position 
of the resonances deduced from the exact solution shown in 
Fig. [3jc). We compare those two quantities for different po¬ 
tential depths (a) and widths (b). 


arguments. For example for a waist of 140 /jl m/32 ~ 4 
fi m, only 3 resonances remains. Remarkably, the local 
envelope approximation at the heart of the semiclassical 
approach still holds when the potential contains only a 
few lattice periods. 


In contrast with the usual optical Fabry-Perot cavity, 
a cavity based on two spatial gaps has a finesse that 
strongly depends on the energy: 


(Je 2 


2 V / 1/T 2 (£) - 1/T(E) 


-1 


(27) 

where A e accounts for the energy difference between two 
successive peaks, and ge is the standard deviation of each 
peak. Figure [5] provides a comparison of the finesse calcu¬ 
lated from Eq. (J2ZD with that obtained from the integra¬ 
tion of the Schrodinger equation. The finesse decreases 
(almost exponentially) with the energy. The quality of 
the agreement is quantitatively evaluated by computing 
the relative error A T{E) between the two approaches. 
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The inset of Fig. [5] shows that this error remains lower 
than 2%, in the considered energy range and tends to 
increase slightly for lower energies. 



E ( E r ) 


FIG. 5: (Color online) Inverse of the finesse as a function 
of the energy, obtained from the exact solution Tpw(F) 
(black dots) and from the semiclassical calculation Tfp(F) 
(green/grey line). Inset: relative error, AF(F) = |Fpw — 
Ffp|/Ffp- The parameters are the same as in Fig. [3] 

The semiclassical expression for the transmission pro¬ 
vides a simple way to evaluate the lifetime of a wave 
packet trapped in the cavity at a mean energy that coin¬ 
cides with that of a resonance. Indeed, the expansion of 
the phase about the resonance (</?(F) ~ i r) reads 


waves. Let us finally emphasize that the level spacing 
and the decay rate can also be engineered for a fixed 
cavity size by renormalization of the mass using an ex¬ 
tra superimposed optical lattice [6]. Such cavities are 
good candidates to observe the atom blockade effect in 
close analogy with its electronic counterpart (Coulomb 
blockade) [33]. 


B. Multiple barrier landscape 




-300 -200 -100 0 100 200 300 


^ (/im) 


FIG. 6: (Color online) Position dependent imaginary part 
of the Mathieu exponent, ft, as a function of the velocity 
(2 E/m ) 1//2 for a lattice with a Gaussian envelope (w — 140 
/im) and for a maximum potential Uo — 10Fr. 


sin (<p(E) + 5<p(E)) « ~g^ SE = —, (28) 

where the semiclassical time r corresponds to the time 
required to travel back and forth in the cavity at that 
energy 


T = -27, P 

Jo 9E 



dz 

v(z;E ) 


(29) 


with v(z;E) = (1 /H)dkE is the semiclassical velocity. 
The expansion of the transmission about the resonance 
reads 


Tpp(F) 


1 + 


4 t 2 
1 - T(F) h 2 



(30) 


We have represented in Fig. [31(b) the position depen¬ 
dent imaginary part of the Mathieu exponent for a lat¬ 
tice having a Gaussian envelope (w = 140 /im) and for a 
maximum potential depth Uo = Fr. By increasing the 
potential depth, we increase the number of lateral cavi¬ 
ties. This is illustrated in Fig. [6l In principle, one could 
implement a more involved envelope to generate a ran¬ 
dom distribution of spatial gap barriers. The physics of 
Anderson localization could therefore be revisited in this 
context jr| . 

However, one has to be careful with the modification of 
the picture for local spatial gaps due to multiple interfer¬ 
ences as already discussed in Sec. IIV Ai To put forward 
this point, we consider the specific case of an optical lat¬ 
tice with a periodic envelope (see Fig. 0(a)): 


from which we infer the time decay, To, 

r °=< si > 

More generally, the size of the cavity provides an en¬ 
ergy scale related to the difference of energy between 
two successive resonances, A#. The finesse of a given 
resonance has an energy width equals to h/ro. For pa¬ 
rameters such that the quality factor, A^ro/S, is large, 
this system realizes a true Fabry-Perot cavity for matter 


V 0 (z) = Uo 


1 + s cos 



(32) 


where Fo>0,F>dis the spatial period of the enve¬ 
lope and e > 0 the amplitude of the spatial modulation 
(see Fig. 0(a) plotted for 5 = 0.25, D = 5 d). A naive 
approach based on the local envelope approach yields a 
periodic pattern for the spatial gaps. Actually, this prob¬ 
lem can be solved exactly by noticing that the expansion 
of the total potential U(z) contains nothing but a super¬ 
position of four periodic potential of spatial periods d, F>, 
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Dd/ (D ±d). For the sake of simplicity, we consider that 
the spatial periods are commensurable. In this case, the 
corresponding band structure can be readily calculated 
using the appropriate Brillouin zone and the polychro¬ 
matic Hill method 0123 • In Fig. [Hb), we have rep¬ 
resented the forbidden band gaps resulting from such a 
calculation as black bands. They are completely flat in 
contrast to the intuition based on the picture proposed 
by the local approximation. This feature results from the 
fact that the envelope is periodic and therefore spatially 
infinite. This shows the limit of the local approach. 



FIG. 7: (Color online) (a) Superlattice potential (l32l) with 
e = 0.25 and D = 5 d. (b) The imaginary part of the wave vec¬ 
tor obtained from the local Mathieu equation is represented 
with the green/grey color code. The exact solution yields flat 
bands represented as black bands. 


C. Phase space engineering 

Spatial gaps enable one to generate non trivial phase 
space correlation. We provide a concrete example based 
on our experimental results. The experiment has been 
described in detail in HU. In brief, a rubidium Bose- 
Einstein condensate is prepared in a magnetic state F = 
1, rriF = —1 in an off-resonance crossed dipole trap. The 
vertical beam of the dipole trap is switched off and the 
BEC is launched in the horizontal optical guide provided 
by the other arm of the crossed dipole trap at a mean ve¬ 
locity v = 8.2 mm/s using a magnetic pulse. A finite size 
optical lattice having a Gaussian envelope (w =140 fi m 
and depth 2.5 Er) and centered about the atomic wave 
packet is ramped on adiabatically in 1 ms just after the 
magnetic launching. The BEC further propagates in the 


lattice while being confined transversally by the horizon¬ 
tal guide. After 20 ms, a wave packet is emitted in the 
opposite direction to that set by the launching proce¬ 
dure. This wave packet originates from the partial tun¬ 
neling through a lateral spatial gap of a wave packet that 
has been previously partially reflected by the other sym¬ 
metric spatial gap (see Fig. Eta)). The size of the cavity 
delimited by the spatial gaps increases with the energy of 
the wave packet. However the corresponding increase in 
velocity is not sufficient to overcome the increase in size 
and the high energy part of the wave packet is emitted 
after the low energy part. As a result the emitted wave 
packet acquires a non trivial phase space correlation (see 
the qualitative picture of this effect in the phase space 
representation of Fig. [8](a)). The experimental signature 
of this correlation is the focussing of the packet at a few 
hundreds of micrometers from the spatial gap followed 
by a subsequent defocussing of the packet. Figure [51(b) 
summarizes our experimental results in which this effect 
can be clearly observed. A good agreement is observed 
with the numerical result that include the finite optical 
resolution of the experiment. The correlation in phase 
space of the emitted packet results from the variation of 
the size of the cavity as a function of the energy, a pa¬ 
rameter that can be tuned experimentally by shaping the 
envelope of the lattice 0. Spatial gaps therefore appear 
here as a way to realize a lens in time, and provide a new 
tool for atom optics. 


V. SPATIAL GAP STRUCTURES IN HIGHER 
DIMENSIONS 


Up to now, we have studied the presence of spatial gaps 
and their applications in one dimension. However, it is 
possible to generalize this concept to higher dimensions. 
In this section, we show how the walls provided by an ap¬ 
propriate shaping of the envelope of an optical lattice can 
be used to design matter wave cavities in higher dimen¬ 
sions, with new dynamical features absent in dimension 
one. 

Consider for instance a 2-dimensional square lattice 
potential generated by two pairs of orthogonal and 
counter propagating laser beams (see Fig. [91(a)). We as¬ 
sume that there is no phase relation between orthogonal 
beams. As a result of the Gaussian envelope of each 
beam, the potential experienced by the atoms reads: 


U(x,y) = -- 


U 0 


■ exp 


v 2 r 


2x 2 




sm 


2 

V d / 


U 0 


- exp 


1 + 


2 y 2 


1 + 


sm' 


■(?) 


(33) 


Figures [91(b) to (e) represent different matter wave cav- 
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2 (jum) 



FIG. 8: (Color online) (a) Sketch of the emission of a wave 
packet through a spatial gap tunnel barrier, (b) Evolution of 
the standard deviation of the size of the emitted wave packet 
as a function of time: Experimental results (black squares), 
numerical result (dashed line). The numerical results take 
into account the finite optical resolution of the experimental 
imaging system (15 /mi). 


ities (generated by the same laser configuration) delim¬ 
ited by spatial gaps for different pseudo-energies and po¬ 
tential depths. To determine the shape of those cavities, 
we proceed in the following manner: (i) at each position 
we calculate the 2D band diagram associated with the 
corresponding local depth, (ii) we fix a pseudo-energy 
Eq , and (iii) we plot a black dot if the pseudo-energy lies 
in a gap. This representation gives a direct insight on the 
shape of the cavities that can be designed in 2D. Actually 
in 2D the Mathieu exponent becomes an anisotropic vec¬ 
tor at each position and cannot therefore be plotted on a 
2D diagram. Interestingly, cavities with different topolo¬ 
gies can be designed. In Fig. [9l we give a few examples 
of simply (d) and non-simply (b), (c) and (e) connected 
cavities depending on the pseudo-energy and potential 
depth. 

Cavities such as the one seen in Fig. [91(d) could be 
used to probe effects of complex dynamics and tunnel¬ 
ing. Indeed, such cavities are akin to billiards with bar¬ 
riers instead of infinite walls. The shape of the cavity is 
known to control the dynamics, allowing the possibility 
to reach chaotic dynamics in parts of the phase space, 


Y(x,y) 



-2 -1.5 -1 -0.5 0 0.5 1 1.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 


x/w x/w 

FIG. 9: (Color online) (a) Sketch of the potential (l33l) for a 
square lattice realized with Gaussian beams (see text). The 
Gaussian parameters are xr/wo — Vr/wo — 20. (b-e) Spatial 
gaps (black) for different values of the pseudo-energy Eo and 
the potential depth Uq\ (b) Eo — 0.05F?r and Uo = Er/2 ( c ) 
Eo = 0.25Er and Uo = Er/2, (d) Eo = 0.35Er and Uo = Er 
(e) Eo = 0.5Er and Uo — Er/2. 


a feature which is not present in dimension one. This 
makes such cavities interesting to simulate the effects of 
chaos with cold atoms in a closed two-dimensional sys¬ 
tem, in contrast with the studies in [37j which concerned 
open systems where the chaotic effects where limited by 
the short time spent in the chaotic region. Beyond the 
study of quantum chaos in a new setting, the presence of 
tunable tunnel barriers at the frontier of the cavity opens 
the way to new types of experimental studies of tunnel¬ 
ing in presence of chaos (chaotic tunneling, chaos assisted 
tunneling) j38|. Such effects have been demonstrated to 
be important in real systems but have been difficult to 
observe experimentally in full detail. The versatility of 
the setup described here could then allow to character¬ 
ize this regime in great detail by varying the different 
parameters of the system. 

We also emphasize that this technique provides a new 
method to generate curved wave guides [39l Idlj as illus¬ 
trated in Fig. [91(d). The geometry of Fig. [91(b) is also 
well adapted for designing matter wave guided structure 
with two symmetric paths for interferometric purposes. 
Interestingly the quality of the symmetry is here auto¬ 
matically ensured by the beam shape. Finally, let us 
point out that the generalization in 3D of cavities delim¬ 
ited by spatial gaps would correspond in the small size 
regime to realizing a quantum dot for atoms, similarly to 
the quantum dots for electrons in mesoscopic physics. 
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VI. TIME-DEPENDENT BAND GAP 
STRUCTURES 

The studies and applications presented in the preced¬ 
ing sections used static potentials. However, it is possible 
to modulate the potentials in time, leading to new effects 
with interesting applications [23|, |42j . In this section, we 
discuss in detail the concept of spatial gaps in an optical 
lattice with a periodic time-dependent amplitude: 

U tm {z-,t ) = [1 + asin 2 (27rz4)] U(z), (34) 

where a < 1 is the strength of the modulation and v 
its frequency. The solution of the Schrodinger equation 
can be expressed locally in terms of the commonly called 
Floquet-Bloch function [23j. Physically, the modulation 
opens new gaps. In the perturbative regime (small am¬ 
plitude of modulation), the modulation at a frequency 
v generates a new gap when hv matches the difference 
A E of energy between two different bands of the band 
diagram obtained in the absence of modulation. We will 
consider the scattering of a non interacting wave packet 
onto the time-modulated potential, studying in turn the 
case of the wave packet initially outside and then inside 
the optical lattice. 

A. Initial wave packet outside the modulated 
optical lattice 

In this subsection, we thus consider a wave packet ini¬ 
tially outside the potential (|M|) with Uo = 2Er and 
a waist w = 140 /im. The incident wave packet has 
a mean velocity v = 0.6ur and a velocity dispersion 
Av = 2.34 x 10 -4 m.s -1 (Av <C v). With such pa¬ 
rameters and in the absence of time modulation (a = 0), 
the wave packet is totally reflected (the probability of 
transmission is lower than 1 %). To further characterize 
the time evolution of the wave packet, it is convenient to 
use the semiclassical picture according to which the wave 
packet is described in terms of a fictitious particle having 
the same mean velocity [23]. 

Far from the potential, the mean energy of the ficti¬ 
tious particle is purely kinetic and equal to mv 2 / 2. This 
sets the value of the pseudo-energy once inside the lattice. 
The adiabatic evolution when the fictitious particle en¬ 
ters the lattice is represented by the dashed black arrows 
(Fig. [TOlf a). (b) and (c)). Choosing an appropriate value 
for the frequency modulation (zq = 0.3 kHz), one can 
promote the particle to the upper band in order to feed 
the main cavity centered about z = 0 (see Fig. [TOl cH. 
Alternatively, one can use a larger frequency (v 2 = 3 
kHz) to drive the fictitious particle into a lower band so 
as to populate a lateral cavity centered about z = —100 
fi m for our parameters (see Fig. ITUlf c)). The resonance 
transitions between bands occur at well-defined positions 
as a result of the envelope and lead to the loading of the 
cavity only in the case of a negative velocity i.e. once the 
atoms have been reflected by the gap since the transition 


connects two energy bands with opposite slopes (the lo¬ 
cal velocity is directly proportional to the local slope of 
the band in the semiclassical picture Q). 

The probability of transition increases with the modu¬ 
lation amplitude a. According to our numerical simula¬ 
tion, for zq = 0.3 kHz and a = 0.5, 10 % of the atoms of 
the incoming wave packet are coupled to the main cav¬ 
ity delimited by the gap between bands II and III. The 
width of the cavity is about 200 fi m for our parameters. 
Figure fTOl d) shows clearly the large oscillation amplitude 
of the packet inside the main cavity. This oscillation is 
accompanied with a focusing of the trajectories that re¬ 
sults from the velocity dependent coupling to the cavity. 
With zq = 3 kHz, one transfers the atoms in the small 
size lateral cavity corresponding to band I (see Fig. Mb) 
and (c)). We observe periodic losses due to the fact that 
the modulation is maintained and drives regularly (at a 
given position) the outcoupling of atoms from the cav¬ 
ity. By stopping the modulation (after 75 ms), one could 
keep many atoms in the lateral cavity for a very long 
amount of time since the spatial gaps are very large and 
the tunneling rate is therefore completely negligible. 


B. Initial wave packet inside the modulated lattice 

We then consider in this subsection the case of a wave 
packet (Fig HTf a)) initially at the center of the time- 
dependent potential dm and with a mean velocity v such 
that the corresponding kinetic energy is equal to the one 
at the middle of the second band in the vanishing po¬ 
tential depth limit: mv 2 /2 = 0.562£r. The potential 
depth is then adiabatically increased up to its final value 
Uo = 2£r (A—>A 5 ), so that the pseudo-energy of the fic¬ 
titious particle is close to -0.3 Er (Fig. [TTI b)). 

In the absence of modulation, the wave packet os¬ 
cillates inside the cavity delimited by the gap between 
bands I and II (Fig.HTlfd)). In the presence of the modu¬ 
lation, the wave packet propagates inside the cavity with 
a fixed pseudo-energy until it reaches the resonant condi¬ 
tion for which it can be promoted to another band. For 
instance, with a modulation frequency v = 5 kHz, the 
transition occurs at B. About 20% of the density goes 
from B to B\ Actually, a more careful analysis reveals 
that two different paths can be followed by the atoms: 
(i) they are promoted on the cavity delimited by the gap 
that separates bands III and IV. Those atoms oscillate 
with the same amplitude but with a smaller time period 
since they have more energy and thus a larger velocity 
than in the original band; (ii) They are transferred to 
band III just below the gap (III to IV) and continue their 
propagation up to point C where a new resonance con¬ 
dition occurs. As a result of this resonance, atoms can 
be transferred (C to C’) to the lateral cavity on the low¬ 
est band I. Part of the atoms continue their propagation 
(C to D) leaving the optical lattice with a mean veloc¬ 
ity 1.2 ur = 8.49 mm/s, larger than v (these processes 
are illustrated in Figs. □He) and UTlf e)). As exemplified 
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FIG. 10: (Color online) (a) Wave packet scattering (mean 
velocity v m 0 . 65 ^r) on an amplitude modulated lattice 
(HU) of depth Uo = 2F?r. (b) Band diagram in the van¬ 

ishing depth limit (green/grey line) and for a local depth 
Vo(z — —104 jLtm) = 0 . 9 F?r (black line). Green/grey dots 
indicate the transfer between bands induced by the propaga¬ 
tion in the presence of modulation. The wave packet mean 
energy of the wave packet is mv 2 /2 = 0 . 42 F?r. The solid 
line black arrows depict the interband resonance transitions. 

(c) Position dependent imaginary part of the Mathieu expo¬ 
nent. The wave-packet scatters on the spatial gap (II —> III), 
(d-e) Propagation of a non-interacting wave packet (mean ve¬ 
locity v , velocity width Av = 2.34 x 10 -4 m/s initially at 
z = —400 jLtm) obtained from a numerical integration of the 
Schrodinger equation based on the split-step Fourier method: 

(d) with a modulation of frequency v\ — 0.3 kHz and (e) 
^2 = 3 kHz (modulation amplitude a = 50 %). 


here, the details of the atomic wave packet evolution can 
be fully understood using the Bloch diagram combined 
with interband transitions instead of the non perturba¬ 
tive Floquet-Bloch diagram approach. 

As illustrated above, the transitions resulting from the 
time modulation of the envelope offer a wide variety of 
possibilities to manipulate propagating wave packets. In 
a scattering experiment with the wave packet initially 
outside the lattice, the output channels can be controlled 
by the modulation. In the case of a wave packet placed 
initially at the center of the optical lattice with a fi¬ 
nite velocity, we have shown how one can load atoms 
in lateral and/or centered cavities delimited by spatial 
gaps. The coupling between the cavities can be engi¬ 
neered by changing the modulation amplitude or/and the 
frequency. Furthermore a wave packet with a non-zero 
momentum trapped in a finite size optical lattice can 
be considered as an atom reservoir. Atoms can then be 


FIG. 11: (Color online) (a) Wave-packet propagation (mean 
velocity v = 0.75i>r) inside an amplitude modulated lattice 
([34]) of depth Uo = 2F?r. (b) ” Position” of the wave packet 
on the band diagram for different potential depths (i.e. dif¬ 
ferent positions z). Adiabatic loading: from Uo = 0 (z = 0) 
(green/grey dashed line) to Uo = 2 E r (black dashed line) in 
t — 1 ms. Propagation: Vo(z cz 40 /im) — Uo cz. 1.7 Er (black 
line) and Vo(z ~ 100 fin i) = Uo — 0.9 Er (green/grey line). 
The dotted arrow represents the evolution of the pseudo¬ 
energy during the adiabatic loading of the wave packet. The 
dashed black arrows represent the propagation of the wave 
packet at a constant pseudo-energy, and the solid black ar¬ 
rows the possible transitions between two bands induced by 
the modulation, (c) Position-dependent imaginary part of 
the Mathieu exponent, (d-e) Propagation of a non inter¬ 
acting wave-packet (mean velocity v, velocity width Av = 
2.34 x 10 4 m/s) initially at z = 0 /im: (d) in the absence 
of modulation and (e) in the presence of a modulation at fre¬ 
quency v — 5 kHz, with a = 33%. 


outcoupled at will by an appropriate choice of the mod¬ 
ulation frequency of the optical lattice depth. Such a 
device realizes a tunable source of atoms with well de¬ 
fined velocity properties, an important prerequisite for 
many quantum devices. 


VII. CONCLUSION AND PERSPECTIVES 

We have studied in detail the concept of spatial gaps 
created by optical lattices with a spatially varying enve¬ 
lope, and presented several applications for matter wave 
engineering and atomtronics. Indeed, we have shown 
that they correspond to tunnel barriers with properties 
difficult to attain by other techniques. In one dimen¬ 
sion, they enable to build Fabry-Perot cavities for matter 
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waves with tunable parameters in direct similarity with 
optics, to design multiple barrier systems and to shape 
the phase space distribution of a wave packet. In higher 
dimensions, we have shown that the generalization of the 
concept of spatial gaps makes possible to engineer curved 
wave guides and cavities similar to quantum dots for mat¬ 
ter waves whose properties can be tuned. Such systems 
correspond to different dynamics as the parameters are 
varied, and enable the study of complex properties absent 
in one dimension such as complex dynamics and chaotic 
tunneling. The time modulation of the optical potential 
adds a new versatility to this kind of systems, enabling 
to load cavities or transfer atoms from one to another 
through the modulation. This could allow a precise en¬ 
gineering of matter wave packets, in particular leading 
to the possibility of new tunable sources of atoms with 
specific velocity properties. 

In addition, spatial gaps provide a test bed for study¬ 
ing the role of dimensionality and interactions on tun¬ 
nel effect [43|, |H|. In particular, the tunnel effect is 
deeply different in higher dimensions, enabling to probe 
the interplay of complex dynamics and tunneling such as 
chaos-assisted tunneling or chaotic tunneling, phenom¬ 
ena which are absent in one dimension. The addition 
of interaction effects has been the subject of very lit¬ 
tle work and could lead to new types of physical effects. 
One could also envision to revisit, in the context of a 
finite-size optical lattice, the dynamical instability that 
have been observed with a Bose-Einstein condensate in 
the Thomas Fermi regime, moving in a lattice [45|. This 
instability can be explained in terms of four wave mix¬ 
ing [46j triggered by the dispersion relation inside the 
lattice [47H49]. As a result of such processes, an atom 
is promoted in the high energy part of the lower band 
while another ends up (by energy conservation) at its 
bottom. Since the tunneling is favored at high energy, 
a cooling by coherent collisional processes could be ob¬ 
served in this manner. Another perspective could be to 
investigate the new possibilities offered by spatial gaps 
using spin-dependent optical lattices [5CH52j. 

We therefore think that the spatial gaps created by 
the local band gap structure of optical lattices enable to 
build new types of physical systems which can serve as 
tools for the developing field of atomtronics and can lead 
to the investigation of new physical effects. 
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Appendix A: Perturbative treatment in the weak 
potential approximation 

Here, we work out a first order perturbation theory 
valid inside and outside the first gap which complements 
Sec. Ill Cl Consider the lowest gap of a periodic potential. 
In the weak potential depth limit, only the first two bands 
matters for the dynamics and the eigenvalue problem is 
equivalent to the two coupled equations [53| : 


I £kVk — ( Uo/4)vk-kB. = E y k 

\-(Uo/4)vk + £k-kn v k-k n = EVk-k n 
The solutions of this system are: 


E\ 


k ^ — 2 fc R ) i 71 


£k £k — k^i 


u 2 

—, (Al) 
16’ V } 


with £k = h 2 k 2 /(2m). Using the notation k = k/k^ — 1/2 
and s = Uo/Er , we can rewrite Eq. m as 


k 

£r 




from which we infer the value of the wave vector as a 
function of the energy 


*'«<*)/**=q ± 

(A2) 

Interestingly this first order approximation formula gives 
the Mathieu exponent in the band gap when an imag¬ 
inary part is present. The relative error e between the 
exact Mathieu exponent and the first order approxima¬ 
tion is defined by 



e 


max 
0 <e<e r 


{ A"exact (E)-K(E) \ 

V ^exact (E) ) • 


(A3) 


The essential contribution to e originates from the en¬ 
ergies close to the gap. We report below (Table [Aj) the 
value of e for different potential depths Uq. 


Uo/E r 

0.1 0.3 0.5 1 

e 

0.008 0.042 0.089 0.276 


Table 1: Numerical values of e for various potential depths 

Uq. 
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Appendix B: Perturbative treatment for the second 

gap 

Here, we extend the treatment of Sec. Una to the sec¬ 
ond gap. The lower border of the second gap is located at 
the center of the Brillouin zone (k = 0). By symmetry, 
we shall evaluate its width by considering three modes 
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vo, v kR and —v kn - The energy at the edge of the gaps 
are obtained at the lowest order by solving the following 
coupled mode equations: 

(■ Er - Uo/2)v-k n - (Uo/tyvo = Ev-k R , 

— (ko/4)r>-/c R — (Uo/2)vq — (Uo/4:)vk R = Ev o. 
~{Uo/4:)vo + (Er - Uo/2)vk R = Evj CR . 

(Bl) 

One finds E_ = Er (1 — 5/2) and E+ = 
Er ^1 — s + \A + s2 /^j /2. In this case, the en¬ 
ergy difference scales as s 2 : (E± — E-)/Er = 
(a/1 + s 2 / 2 — l)/2 ~ s 2 /8. To determine the imaginary 
part of the wave vector inside the second gap, we shall 
follow the same method as for the first gap (see Sec. hej. 
For a given energy E x = (E+ + E-)/ 2 + #(.E + — EJ)/2 
with — 1 < x < 1 inside the gap, we search for a 
solution involving the three following wave vectors 
k± = ±(1 ± iK/2)k R and go = ik R K/2: 

((1 - iK/2) 2 - s/2 - e x )v- kR - (s/4)v 0 = 0, 
~(s/4)v-k R + {-K 2 / 4 - s/2 - e x )v 0 - (s/4)v fea = 0, 
— (s/4)v 0 + ((1 + iK/2) 2 - s/2 - e x )v kR = 0, 

(B2) 

with e x = E x / Er . The determinant of this three-mode 
system yields a polynomial in K of order 6 whose unique 


positive real solution is the imaginary part we are inter¬ 
ested in. In Fig.[l2j we represent the corresponding result 
for s = 0.5. We have also plotted the result of a five-mode 
approximation which is in very good agreement with the 
exact result. 



E/Er 

FIG. 12: Imaginary part of the wave vector normalized to 
/cr/2, K(E,s), in the second band gap for s = 0.5: exact 
calculation (solid line), perturbative calculation to the lowest 
order (dashed line) and perturbative calculation pushed to 
the next order (black square). 
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